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Section  1 
Introduction 


This  report  covers  two  major  topics: 

1.  modeling  of  infrasound  signals  from  atmospheric  explosions  and  evaluation  of  International 
Monitoring  System  network  performance  using  data  from  historic  Soviet  nuclear  tests  and 
other  infrasound  data  sources;  and 

2.  analysis  of  infrasound  instrumentation  through  a  program  of  experimentation  and  theoretical 
modeling. 

The  Institute  for  the  Dynamics  of  the  Geospheres  (IDG)  in  Moscow,  Russia,  has  an  archive  of 
infrasound  recordings  from  Soviet  atmospheric  nuclear  tests  that  were  conducted  in  1957  and 
1961.  IDG  has  digitized  a  total  of  220  infrasound  waveforms  from  22  atmospheric  nuclear  tests. 
Two  events  were  high  altitude  explosions  at  Kapoustin  Yar,  8  others  were  at  the  Shagan  River 
Test  Site,  and  the  other  12  were  at  Novaya  Zemlya.  138  of  the  waveforms  have  measurable, 
unclipped  signals,  known  instrument  responses,  yields,  and  calibrations.  We  have  been  modeling 
this  data  in  two  ways:  first,  by  using  the  data  to  place  constraints  on  infrasound  scaling  relations; 
and  second,  by  numerical  modeling  of  the  infrasound  signals.  Scaling  relations  estimate  the 
pressure  as  a  function  of  yield  and  range.  Several  scaling  relations  exist  in  the  literature,  and  the 
differences  between  them  lead  to  large  differences  in  estimates  of  the  detection  threshold  of  the 
International  Monitoring  System  (IMS).  We  find  that  the  data  is  in  best  agreement  with  a  scaling 
and  attenuation  relation  developed  by  Los  Alamos  National  Laboratory  (LANL)  which  can  be 
written  as  log  P  =  3.37  +  0.68  log  W  - 1.36  log  R  where  P  is  zero  to  peak  pressure  amplitude  in 
Pascals,  W  is  the  yield  in  kilotons,  and  R  is  the  source  to  receiver  distance  in  kilometers. 

IMS  network  detection  capability  is  predicted  using  the  network  simulation  programs  NetSim 
and  XNICE.  NetSim  uses  joint  probabilities  together  with  a  scaling  relation  to  calculate  the 
detection  threshold.  XNICE  has  the  additional  capability  of  estimating  the  detection  threshold  for 
non-Gaussian  noise  distributions.  We  find  that  IMS  noise  levels  have  a  time-dependent  and  non- 
Gaussian  distribution.  The  noise  is  characterized  by  a  fairly  well-defined  minimum  level,  but 
highly  variable  maxima.  The  two-station  90%  probability  IMS  detection  threshold  is  found  to  be 
about  one  kiloton  with  higher  thresholds  in  ocean  regions  and  lower  thresholds  on  land. 

A  typical  infrasound  recording  instrument  uses  an  array  of  pipes  connected  to  a  central  manifold 
to  spatially  filter  random  pressure  fluctuations  along  the  ground  and  so  enhance  the  signal-to- 
noise  ratio  (S/N).  We  present  a  rapid,  accurate,  method  of  estimating  the  performance  of  any  pipe 
configuration  used  for  infrasound  recording.  We  first  describe  experimental  work  that  provides 
the  physical  basis  for  simulations  and  permits  estimation  of  parameter  values.  We  then  describe 
the  propagator  method  used  for  the  calculations.  The  performance  of  a  pipe  array  depends  on  the 
character  of  the  noise,  and  so  meaningful  modeling  of  infrasound  pipe  array  performance  hinges 
on  accurate  knowledge  of  the  temporal  and  spatial  distribution  of  pressure  fluctuations  over  the 
array.  We  discuss  the  physical  basis  for  the  noise  model,  and  limitations  imposed  by  necessary 
simplifications.  We  present  the  amplitude  and  phase  responses,  and  their  sensitivities  to  different 
parameters,  of  two  different  types  of  pipes  used  to  construct  the  arrays:  closed-ended  pipes  with 
acoustic  inlets  along  their  lengths  versus  open  ended,  otherwise  impermeable  pipes.  Finally,  we 
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simulate  the  S/N  enhancement  for  six  different  pipe  configurations  that  have  been  considered  for 
the  IMS  infrasound  stations,  and  discuss  the  results'  implications  for  instrument  design. 

We  find  that  open-ended  pipes  are  preferred  over  closed-ended  pipes  with  high  impedance 
acoustic  inlets.  The  open-ended  pipes  have  essentially  flat  amplitude  responses,  while  the  closed- 
ended  pipes  act  as  lowpass  filters.  The  open-ended  pipes  also  have  flat  phase  responses,  while  the 
phase  responses  of  the  closed-ended  pipes  vary  much  more  across  the  spectrum  of  interest.  The 
most  important  design  criterion  for  enhancing  S/N  is  having  sufficient  spatial  sampling.  Extra 
ports,  however,  provide  no  advantage  if  spaced  more  closely  than  the  noise  correlation  length. 

All  configurations  performed  well,  and  similarly,  at  high  frequencies,  as  the  high  frequency  noise 
has  the  least  spatial  correlation.  The  differences  in  performance  are  very  significant  at  the  lower 
frequencies.  The  two  70  m  diameter  configurations  provide  the  best  performance  in  our 
simulations  of  S/N  enhancement,  and  their  performances  were  nearly  identical  even  though  one 
had  144  ports  and  the  other  had  only  96.  The  18  m  diameter  configuration  performed  the  worst  at 
low  frequency.  The  most  important  step  that  can  be  taken  now  to  improve  instrument 
performance  is  to  obtain,  and  utilize  in  such  modeling,  more  complete  and  accurate  noise 
models. 
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Section  2 

Infrasound  Scaling  and  Attenuation  Relations  and  IMS  Detection  Capability 

The  International  Monitoring  System  (IMS)  specified  by  the  Comprehensive  Nuclear-Test-Ban 
Treaty  (CTBT)  will  include  sixty  infrasound  stations  designed  to  detect  atmospheric  nuclear 
explosions  (Figure  1).  A  design  goal  is  that  the  system  be  able  to  detect  and  locate  explosions  as 
small  as  one  kiloton  anywhere  in  the  world.  In  order  to  estimate  the  capability  of  the  infrasound 
network,  it  is  necessary  to  be  able  to  predict  the  amplitude  of  an  infrasound  signal  at  any 
location,  and  to  evaluate  whether  the  signal  would  be  detectable  above  noise  levels  at  the 
recording  stations.  Scaling  and  attenuation  relations  are  empirical  and/or  theoretical  equations 
that  relate  the  amplitude  and  period  of  infrasound  signals  to  the  explosion  yield  and  source  to 
receiver  distance.  Several  different  relations  have  been  developed  based  on  theoretical  infrasound 
modeling,  and  on  recordings  of  atmospheric  nuclear  and  chemical  explosions.  The  different 
scaling  relations,  however,  imply  very  different  detection  threshold  levels.  In  the  following 
section  we  review  the  different  relations  and  their  implications  for  detection  thresholds. 


Figure  1.  Stations  of  the  planned  IMS  infrasound  network. 


2.1  Infrasound  Scaling  and  Attenuation  Relations. 

Pierce  and  Posey  (1971)  developed  a  solution  for  the  excitation  of  the  Lamb  edge  mode  as  an 
approximation  for  infrasound  signals  generated  by  atmospheric  tests.  They  showed  that  this 
relatively  simple  approximation  gave  waveforms  that  agreed  well  with  the  first  few  cycles  of 
observed  waveforms  from  multi-megaton  nuclear  tests.  From  this  they  developed  the  relation 
between  the  yield  W  in  kilotons,  zero  to  peak  pressure  P  in  Pascals,  the  period  of  the  first  cycle 
of  the  waveform  T,  and  the  angular  distance  A  given  by: 

IV  =  0.494  P(sin  A)  2  Hs  (cT)V2  ( 1 ) 
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where  c  is  the  sound  speed,  Hs  is  the  atmospheric  scale  height  c1  /yg,  where  yis  the  adiabatic 
expansion  constant,  approximately  1.4  for  air,  and  g  is  the  acceleration  of  gravity.  Posey  and 
Pierce  (1971)  showed  that  this  relation  agreed  well  with  a  data  set  of  observations  of  pressure 
and  period  measurements  from  large  atmospheric  nuclear  explosions. 

The  derivation  of  Equation  1  uses  the  approximation  that  the  period  being  measured  is  much 
longer  than  the  source  duration.  Pierce  and  Posey  estimated  the  source  duration  to  be 
approximately  Ts  =  0.33WI/3  seconds,  which  is  only  about  13  seconds  even  for  a  58  megaton 
explosion,  and  therefore  always  much  less  than  the  measured  infrasound  period.  This 
approximation  has  the  effect  of  making  period  independent  of  yield  and  proportional  to  the  cube 
root  of  distance.  For  the  typical  sound  velocity  (.318  km/sec)  and  scale  height  (8  km)  referenced 
in  the  paper,  the  period  of  the  first  cycle  of  the  Lamb  wave  can  be  written  as 

T  =  l3ARm  (2) 

where  R  is  the  source  to  receiver  distance  in  kilometers  and  T  is  in  seconds.  Equation  2  can  then 
be  combined  with  equation  1  to  give 

W  =  34.8PR1/2(sin  A)'/2  (3) 

The  pressure  is  therefore  predicted  to  be  proportional  to  yield  at  a  fixed  distance  and  to  decrease 
approximately  inversely  with  distance.  As  can  be  seen  from  the  data  in  the  following  section, 
these  predictions  are  inconsistent  with  the  data  which  scales  much  more  slowly  than  linearly  with 
yield  and  exhibits  a  clear  frequency  dependence  with  yield. 

Equation  1  has  been  widely  used  (McKisic,  1997)  and  turns  out  to  be  in  good  agreement  with  the 
data  when  the  observed  period  is  used,  even  though  the  usage  is  inconsistent  with  the  derivation 
of  equation  1.  Blandford  and  Clauter  (1995)  modified  equation  1  by  assuming  that  the  period  T  is 
proportional  to  the  cube  root  of  yield  and  independent  of  distance.  From  this  it  follows  that 
pressure  is  proportional  to  the  square  root  of  yield  and  decays  as  the  square  root  of  distance  from 
the  source.  Blandford  and  Clauter  (1995)  used  the  square  root  of  distance  decay  at  distances  less 
than  20  degrees,  but  modified  the  attenuation  relation  to  have  an  exponential  form  at  distances 
greater  than  20  degrees  to  match  the  data  set  of  Wexler  and  Hass  (1962).  Using  a  larger  data  set, 
keeping  the  same  yield  scaling,  but  modifying  the  attenuation  relation,  Clauter  and  Blandford 
(1998)  derived  the  relation: 

log  P  =  0.92  +  05  log  W  - 1.47  log  A  (4) 

where  A  is  the  source  to  receiver  distance  in  degrees,  W  is  the  yield  in  kilotons,  and  P  is  the  zero 
to  peak  pressure  amplitude  in  Pascals.  They  showed  that  this  relation  was  consistent  with  a 
historical  data  set  of  nuclear  and  chemical  explosions. 

Whitaker  (1995)  derived  the  following  relation  based  on  wind-corrected  infrasound 
measurements  from  Los  Alamos  chemical  explosions: 

Pc=235xlO\R/W',2yU6  (5) 


! 
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where  Pc  is  the  wind  corrected  zero  to  peak  pressure  in  Pascals,  R  is  the  source  to  receiver 
distance  in  kilometers,  and  W  is  the  yield  in  kilotons  (the  original  reference  stated  that  yield  is  in 
tons,  but  should  be  kilotons  (Whitaker,  personal  communication)).  The  wind  correction  applied 
in  equation  5  is 

p=l0-OOl9Vp  (6) 

where  v  is  the  component  of  the  stratospheric  wind  velocity  in  meters/second  in  the  direction  of 
wave  motion  and  P  is  the  measured  pressure. 

Russian  scientists  at  IDG  have  used  the  following  relationships  to  make  approximate  yield 
estimates: 


Pm  =  kpWm!R 

(7) 

T  =  kTWV3  /  R 

(8) 

(9) 

where  Pm  is  the  zero  to  peak  pressure  in  Pascals,  T  is  the  signal  duration  from  the  first  arrival  to 
the  moment  when  the  signal  degrades  to  noise  level,  and  90%  of  the  signal  energy  is  contained  in 
the  frequency  band  below  angular  frequency  0)m.  The  three  constants  are  empirically  determined 
and  are  given  in  Table  1. 


Table  1.  Russian  scaling  constants 


kp  Pa-km/kt1/3 

kT  s/ktl/3km 

km  kt1/3/s 

Downwind 

2000 

0.050 

1.6 

Perpendicular  to  wind 

1000 

0.025 

1.2 

The  equations  above  can  be  written  in  a  consistent  format  as: 


log  P  =  -L54  +  log  W  -  0.5 log(R  sin  A) 

Pierce  and  Posey  (1971) 

(10) 

log  P  =  0.92  +  0.5  log  W  - 1.47  log  A 

AFTAC  (Clauter  and  Blandford,  1998) 

(ID 

log  P  =  3.37  +  0.68  log  W  - 1.36  log  R 

LANL  (Whitaker,  1995) 

(12) 

log  P  =  3.00  +  0.33  log  W  -  log  R 

Russian  -  Crosswind 

(13) 

log  P  =  3.30  +  0.33  log  W  -  log  R 

Russian  -  Downwind 

(14) 

where  P  is  zero  to  peak  pressure  in  Pascals,  W  is  yield  in  kilotons,  R  is  distance  in  kilometers, 
and  A  is  distance  in  degrees. 

These  scaling  relations  have  very  different  implications  for  the  infrasound  detection  threshold. 
Yield  estimates  for  the  threshold  pressure  level  differ  by  several  orders  of  magnitude,  even 
though  each  relation  was  constrained  by  some  infrasound  data  set.  Table  2  shows  the  calculated 
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yield  at  a  nominal  detection  threshold  of  0. 1  Pascal  for  each  of  the  scaling  relations.  The  most 
important  factor  is  the  exponent  in  the  pressure/yield  relation.  With  the  Pierce/Posey  relation, 
which  has  a  yield  exponent  of  1,  pressure  drops  off  much  more  rapidly  with  yield  than  with  the 
other  relations  leading  to  a  very  high  threshold  level.  The  0.33  yield  exponent  in  the  Russian 
relations,  however,  implies  a  very  slow  decrease  in  pressure  with  yield  and  leads  to  very  low 
threshold  levels.  The  Whitaker  and  Clauter/Blandford  relations,  which  have  yield  exponents  of 
0.68  and  0.5,  respectively,  predict  intermediate  threshold  levels. 

Table  2.  Detection  capability  for  a  nominal  detection  threshold  of  0.1  Pascal.  The  table  shows  yield  in  kilotons  for 
each  scaling  relation  and  source  to  receiver  distances  from  10-60  degrees. 


10 

20 

30 

40 

50 

60 

Clauter/Blandford 

0.13 

0.97 

3.2 

7.4 

14 

24 

Whitaker 

0.46 

1.8 

4.2 

7.4 

11 

17 

Russia  Crosswind 

0.0014 

0.01 

0.037 

0.088 

0.17 

0.30 

Russia  Downwind 

0.0002 

0.0014 

0.0046 

0.011 

0.02 

0.037 

Pierce/Posey 

47 

93 

140 

180 

220 

260 

2.2  Soviet  Infrasound  Data. 

IDG  has  an  archive  of  approximately  300  recordings  from  34  Soviet  atmospheric  nuclear  tests 
that  were  conducted  in  1957  and  1961.  20  of  these  explosions  were  located  at  the  Novaya  Zemlya 
test  site,  12  at  Semipalatinsk,  and  2  at  Kapoustin  Yar  (see  Figure  2).  Of  these,  220  recordings 
from  22  of  the  tests  recorded  at  stations  from  1000  to  5000  km  were  found  to  be  of  adequate 
quality  for  digitization.  The  yields  of  these  tests  range  from  0.4  KT  to  58  MT.  The  data  set 
includes  two  high  altitude  explosions  and  the  largest  (58  megaton)  atmospheric  explosion  ever 
detonated.  The  explosions  corresponding  to  this  data  set  are  listed  in  Table  3.  The  number  of 
records  listed  in  the  table  is  the  total  number  of  records  digitized  for  each  event. 


20*  30*  40*  50*  60*  70*  80*  90*  100*  11 0*1 20*  130*  140*  150*  160*  170* 


Figure  2.  Stations  recording  inffasound  signals  from  atmospheric  explosions  at  three  Soviet  nuclear  test  sites 
Semipalatinsk  (STS),  Kapoustin  Yar  (KY),  and  Novaya  Zemlya  (NZ). 
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Table  3.  Soviet  atmospheric  nuclear  explosions  and  the  number  of  records  for  each  event  that  have  been  digitized 
by  IDG.  STS  is  the  Semipalatinsk  test  site,  KY  is  Kapoustin  Yar,  and  NZ  is  Novaya  Zemlya. 


Test  Site 

Explosion 

Number 

Date 

Time 

(Moscow) 

Latitude 

Longitude 

Height  of 
Burst  (m) 

Yield 

(KT) 

Number 

Records 

STS 

088 

1961/09/04 

08:00:27 

50.45 

77.74 

725 

9 

2 

STS 

089 

1961/09/05 

09:00:05 

50.45 

77.74 

500 

16 

3 

KY 

091 

1961/09/06 

48.45 

44.30 

22700 

11 

6 

NZ 

095 

1961/09/10 

12:00:14 

73.52 

54.30 

2000 

2700 

23 

NZ 

099 

1961/09/12 

13:08:00 

73.52 

54.30 

1190 

1150 

23 

NZ 

102 

1961/09/14 

12:56:16 

73.52 

54.3 

1700 

1200 

17 

STS 

103 

1961/09/14 

08:59:59.4 

50.35 

77.82 

0.5 

0.4 

1 

NZ 

111 

1961/09/20 

11:12:12 

73.52 

54.30 

1600 

450-20001 

6 

STS 

112 

1961/09/21 

17:01:01.6 

50.33 

77.70 

110 

0.8 

2 

NZ 

113 

1961/09/22 

11:11:00 

73.52 

54.30 

1300 

260 

4 

STS 

114 

1961/09/26 

10:01:19.8 

50.45 

77.75 

665 

1.2 

2 

NZ 

116 

1961/10/02 

13:30:50 

73.92 

54.55 

1500 

250 

13 

STS 

117 

1961/10/04 

10:01:19.9 

50.44 

77.76 

605 

13 

1 

KY 

119 

1961/10/06 

48.45 

44.30 

41300 

40 

4 

NZ 

120 

1961/10/06 

10:00:08 

73.52 

54.30 

2700 

4000 

20 

STS 

123 

1961/10/17 

10:00:00.8 

50.45 

77.75 

505 

6.6 

2 

STS 

124 

1961/10/19 

08:30:42.6 

50.45 

77.73 

710 

10 

2 

NZ 

125 

1961/10/20 

11:07:03 

73.52 

54.30 

10002 

1450 

17 

NZ 

126 

1961/10/23 

11:31:22 

73.5 

54.3 

3500 

12500 

21 

NZ 

128 

1961/10/25 

11:31:05 

73.52 

54.3 

1450 

300 

10 

NZ 

133 

1961/10/30 

11:33:27 

73.52 

54.30 

4000 

58000 

21 

NZ 

147 

1961/11/04 

10:20:23 

73.5 

54.3 

1750 

150-1500 

20 

Seventeen  stations  recorded  data  from  these  tests.  The  stations  are  listed  in  Table  4  and  Figure  2 
shows  the  paths  from  each  event  to  each  station  in  the  data  set.  Absolute  times  are  not  known  for 
any  of  the  waveforms.  Data  were  recorded  on  instruments  with  varying  low  and  high  pass  filters. 
Figure  3  shows  the  response  of  the  three  major  instrument  types.  EDG  has  put  considerable  effort 
into  identifying  the  instrument  parameters  that  were  in  use  at  each  of  the  recording  stations.  The 
instrument  parameters  are  defined  and  described  for  each  station  in  Appendix  A.  Appendix  B 
lists  the  instrument  parameters  for  each  waveform  together  with  IDG  measurements  of  a  number 
of  parameters. 


1  Yield  information  varies.  2000  kt  was  used  in  the  analysis  in  this  report.  IDG  used  450  kt  in  Appendix  B. 

2  IDG  estimate. 
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/Czcn/iumybuo  -  vavmomme  zapartmepucmurfu: 


Figure  3.  Three  major  instrument  types  used  for  recording  infrasound  data. 


Table  4.  Stations  recording  Soviet  atmospheric  tests. 


Station 

Num. 

Station  Name 

Latitude 

Longitude 

Novaya  Zemlya 
distance  (km) 

Semipalatinsk 
distance  (km) 

Kapoustin  Yar 
distance  (km) 

1 

Dubna 

56.7 

37.3 

2700 

900 

2 

Kirov 

58.6 

49.7 

1700 

2100 

1100 

3 

Esheri 

43.1 

40.9 

3500 

2900 

800 

4 

Podolsk 

55.5 

37.3 

2200 

2700 

900 

5 

Oussouriysk 

43.9 

132.0 

5000 

4100 

6300 

6 

Petropavlovsk 

53.1 

158.8 

4850 

5300 

7100 

7 

Yuzno-Sakhalinsk 

47.0 

142.8 

5000 

4600 

6700 

8 

Zvanets 

48.5 

26.4 

3150 

3650 

1400 

9 

Mayly-say 

41.3 

72.5 

3800 

1140 

2250 

10 

Semipalatinsk 

50.4 

80.2 

2850 

- 

2500 

11 

Baykal 

51.8 

107.6 

3450 

2000 

4250 

12 

Olenegorsk 

68.1 

33.3 

1000 

3100 

2300 

13 

Kooldour 

49.2 

131.8 

4450 

3800 

5900 

14 

Kazakhstan 

47.5 

63.5 

2950 

1100 

1250 

15 

Kap  Yar 

48.7 

46.4 

2800 

2250 

- 

16 

Yeniseysk 

58.5 

92.2 

2350 

1300 

3300 

17 

Leningrad  Reg. 

59.9 

30.2 

1850 

2800 

1650 

Some  examples  of  infrasound  records  are  shown  in  Figure  4.  These  records  are  from  tests  that 
range  from  8  kilotons  to  58  megatons  in  yield.  All  of  these  traces  were  recorded  on  instrument 
#2  except  those  for  event  133,  which  were  recorded  on  instrument  #1.  The  frequency  content  of 
the  signals  changes  dramatically  over  this  range. 
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9  KT 
1140  km 


11  KT 

High  Altitude 
1100  km 

250  KT 
4850  km 

2.7  MT 
4850  km 

2.7  MT 
5000  km 


58  MT 
1700  km 

58  MT 
5000  km 


Figure  4.  Infrasound  signals  from  Soviet  atmospheric  tests.  The  top  record  is  from  test  088  at  station  9,  the  second 
is  from  the  high  altitude  test  091  at  station  2,  the  third  is  from  the  test  1 16  at  station  6,  the  fourth  and  fifth 
are  from  test  095  at  stations  6  and  7,  and  the  last  two  records  are  from  the  test  133  at  stations  2  and  7. 
Amplitudes  are  in  Pascals. 


2.3  Measurement  of  Russian  Data. 

All  of  the  Russian  data  were  carefully  measured  in  a  consistent  manner.  The  data  were  first 
filtered  to  remove  long  and  short  period  noise  outside  the  frequency  band  of  the  data.  A 
Butterworth  filter  was  used  with  comer  frequencies  of  .01  and  .2  Hz  for  events  with  yield  less 
than  100  kilotons,  .002  and  .1  Hz  for  events  with  yields  between  100  kilotons  and  2  megatons, 
and  .001  and  .1  Hz  for  events  with  yield  greater  than  2  megatons.  The  amplitude  and  period  were 
measured  as  half  the  maximum  peak  to  peak  amplitude  and  twice  the  time  difference  between  the 
peak  and  trough,  respectively.  Measurements  were  made  on  both  the  acoustic  wave  and  the  low 
frequency  Lamb  wave  if  possible.  Only  data  with  known  instrument  responses  were  measured, 
and  a  digital  correction  for  the  instrument  response  was  made  at  the  observed  period.  A  few 
signals  with  apparent  calibration  errors  were  not  used.  The  final  result  is  that  measurements  were 
made  on  a  total  of  133  waveforms.  107  acoustic  waves  from  17  events  and  96  Lamb  waves  from 
9  events  were  measured. 
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Infrasound  Measurements  Infrasound  Measurements 


10  100  1000  10  100  1000  100C 

Scaled  Range  (R/W*0 .5)  Scaled  Range  (R/W'O  34) 


Infrasound  Measurements  Infrasound  Measurements 


1  10  100  10i  1UU  lmJ 

Scaled  Range  (RAW)  Scated  R8nB® 


Figure  5.  Comparison  of  Russian  data  with  scaling  relations.  Top  left  is  the  LANL  relation  (Whitaker)  Equation  12, 
top  right  the  AFT  AC  relation  (Clauter  and  Blandford)  Equation  11,  bottom  left  is  the  Pierce/Posey 
relation  Equation  10,  and  bottom  right  are  the  Russian  relations  Equations  13  and  14. 

Figure  5  shows  a  comparison  between  pressure  measurements  made  from  the  Russian  data  and 
Equations  10-14.  Also  shown  on  the  figure  are  measurements  from  the  US  Project  Dominic  tests. 
It  is  not  possible  to  put  them  all  on  the  same  plot  because  the  scaling  relations  have  different 
functional  forms.  All  pressures  are  zero  to  peak  amplitudes  in  Pascals.  The  LANL  relation 
(Equation  12)  appears  to  fit  the  data  very  well  over  the  entire  scaled  range,  although  there  is 
considerable  scatter  about  the  line.  The  pressure  measurements  have  not  been  wind  corrected. 

The  AFT AC  relation  (Equation  1 1)  fits  the  lower  yield  data  quite  well.  The  Pierce/Posey  relation 
(Equation  10)  is  a  fairly  good  fit  to  the  high  yield  Lamb  wave  data,  but  does  not  fit  the  acoustic 
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wave  data,  particularly  for  the  lower  yield  events.  The  Russian  relations  (Equations  13-14)  also 
fit  the  data  fairly  well  with  the  crosswind  equations  matching  the  lower  amplitude  data  and  the 
downwind  equations  matching  the  higher  amplitude  data,  however  the  observed  data  falls  below 
the  predicted  curves  for  larger  scaled  ranges,  and  the  data  points  for  the  largest  yield  events  are 
well  above  the  curve,  while  the  lower  yield  events  lie  below  the  curve.  We  conclude  from  this 
that  the  pressure/yield  slope  of  0.33  used  in  the  Russian  relation  is  too  small,  and  the 
pressure/yield  slope  of  1.0  used  in  the  Pierce/Posey  relation  is  too  large.  The  data  is  consistent 
with  the  intermediate  slopes  of  0.5-0.68  of  the  AFTAC  and  LANL  relations,  with  the  LANL 
relation  fitting  the  data  over  the  widest  range. 


Infrasound  Measurements 


Figure  6.  Period  vs.  explosion  yield  for  all  data. 

Figure  6  shows  the  measured  period  plotted  vs.  yield  for  the  Russian  data,  together  with  lines 
showing  cube  root  scaling  of  period  with  yield  with  the  best  fit  to  the  acoustic  and  Lamb  data 
taken  as  independent  data  sets.  The  cube  root  scaling  approximately  fits  the  data,  however  there 
is  a  large  amount  of  scatter  and  the  slow  increase  of  period  with  yield  would  cause  considerable 
uncertainty  in  a  yield  estimate  based  on  measured  period. 

2.4  Infrasound  Magnitudes. 

In  the  analysis  above,  we  used  a  data  set  of  infrasound  waveforms  from  Soviet  atmospheric  tests 
ranging  in  yield  from  6  kilotons  to  58  megatons  to  place  constraints  on  infrasound  scaling 
relations  and  to  estimate  the  detection  threshold  of  the  future  International  Monitoring  System. 
Analysis  of  133  waveforms  shows  that  measured  pressures  are  consistent  with  yield  and 
attenuation  scaling  relations  developed  at  LANL  for  HE  tests,  and  also  fairly  consistent  with  a 
scaling  relation  developed  by  AFTAC.  Because  the  LANL  relation  is  consistent  with  data  over  a 
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very  wide  yield  range,  it  has  recently  been  adopted  as  the  basis  for  an  infrasound  magnitude  by 
the  International  Data  Center  (EDC)  (Brown,  1999).  A  magnitude  is  a  useful  quantity  for  giving 
an  estimate  of  source  size  that  is  independent  of  the  distance  at  which  the  signal  is  measured.  The 
magnitude  equation  is: 


M.  -  log10  P  +  1.361og10  R  -  0.019v  (15) 

where  the  last  term  corrects  for  wind  as  discussed  earlier.  Infrasound  magnitudes  for  the  Russian 
data  set  (without  wind  correction),  are  listed  in  Table  5.  Magnitudes  for  the  acoustic  wave  and 
Lamb  wave  are  listed  separately.  Figure  7  shows  Mj  plotted  vs.  Yield  for  the  seven  Soviet 
explosions.  Also  shown  is  the  LANL  relation,  Equation  12,  rewritten  as  a  magnitude/yield 
relation: 


M,  =  0.68  log  W  +  3.37 


(16) 


As  can  be  seen  in  Figure  7,  Equation  16  fits  the  data  over  this  very  wide  yield  range  quite  well. 


Table  5.  Infrasound  magnitude  for  eighteen  Soviet  explosions. 
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Figure  7.  Infrasound  magnitude  plotted  vs.  yield  for  eighteen  Soviet  explosions. 

2.5  Spectral  Measurements. 


Some  additional  insight  into  the  scaling  laws  can  be  obtained  by  examining  the  spectra  of  arrivals 
with  different  yields  at  the  same  station.  Station  7,  at  Yuzno-Sakhalinsk,  recorded  four 
atmospheric  explosions  with  yields  of  250,  1450,  2700,  and  58000  kilotons.  The  spectra  of  these 
four  arrivals  are  shown  in  Figure  8.  The  shape  of  the  spectra  change  dramatically  over  this  yield 
range,  with  much  more  low  frequency  energy  at  higher  yields. 


Infrasound  Spectra  at  Station  7 


Frequency  (Hz) 


Figure  8.  Spectra  of  four  infrasound  signals  recorded  at  station  7  with  yields  of  250,  1450, 2700,  and 
58000  kilotons. 
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Figure  9  shows  the  change  in  amplitude  for  three  frequencies  plotted  as  a  function  of  yield, 
showing  that  the  slope  of  the  log  amplitude  vs.  log  yield  curve  changes  as  a  function  of 
frequency.  Figure  10  shows  the  slope  of  the  amplitude  vs.  yield  curve  (as  shown  in  Figure  9) 
plotted  as  a  function  of  frequency.  This  figure  shows  that  the  slope  of  the  amplitude/yield  curve 
is  strongly  frequency  dependent,  and  that  the  slope  is  close  to  1,  as  in  the  Pierce/Posey  model,  at 
very  low  frequencies,  but  declines  to  approximately  1/3,  as  in  the  Russian  scaling  laws,  at  higher 
frequencies. 


Log  Yield  (KT) 

Figure  9.  Spectral  amplitude  plotted  vs.  yield  at  frequencies  of  0.005  Hz  (top),  0.015  Hz  (middle),  and  0.050  Hz 
(bottom).  The  slope  of  the  amplitude  yield  curve  decreases  with  increasing  frequency. 


14 


Slope  of  Amplitude/Yield  Curve 


Figure  10.  The  slope  of  the  log  amplitude  vs.  log  yield  curve  plotted  as  a  function  of  frequency. 

The  strong  dependence  of  the  infrasound  spectra  on  yield  suggests  that  the  explosion  source 
function  is  much  longer  than  the  duration  of  the  explosion  and  its  associated  near  field  nonlinear 
effects.  As  noted  earlier,  these  effects  last  only  a  few  seconds,  while  the  spectral  shape  is  affected 
at  periods  of  hundreds  of  seconds. 

2.6  Network  Detection  Simulations. 

In  order  to  predict  the  performance  of  the  proposed  60  station  IMS  infrasound  network,  we 
modified  the  network  simulation  program  NetSim  (Sereno  et  al.,  1990)  to  include  the  models  of 
infrasound  propagation  described  in  Equations  11  and  12.  NetSim  uses  these  equations  to 
calculate  the  pressures  as  a  function  of  yield  and  range.  These  are  used  together  with  station 
locations,  a  noise  model,  a  minimum  signal  to  noise  ratio  for  reliable  measurement,  and  the 
number  of  stations  required  for  a  signal  to  be  reported,  to  determine  the  network  detection 
threshold  as  a  function  of  position  on  the  earth.  We  have  calculated  detection  thresholds  for  the 
proposed  IMS  network  using  station  dependent  noise  level  estimates.  The  calculations  were 
performed  using  the  following  parameters: 

1.  We  used  noise  estimates  from  Blandford,  et  al  (1995),  which  are  based  on  wind 
measurements.  The  log  noise  levels  (Figure  11)  varied  from  -1.37  to  -0.09  (Pascals),  with 
the  highest  levels  being  in  oceanic  regions.  Log  standard  deviation  was  taken  to  be  0.37. 

2.  Simulations  were  performed  for  two  minimum  signal  to  noise  ratio  levels:  2.0  and  1.5. 

3.  Two  stations  detect  infrasound  signals  at  a  90%  confidence  level. 

4.  Four  element  infrasound  arrays  increase  signal  to  noise  ratio  by  a  factor  of  2. 

5.  Propagation  error  has  a  log  standard  deviation  in  log  signal  of  0.3. 

6.  Station  reliability  is  95%. 
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k>gi O(noise),  Pascals 

Figure  11.  Frequency  distribution  of  noise  for  the  IMS  network  estimated  from  wind  measurements  by  Blandford  et 
al.  (1995).  Vertical  axis  is  the  number  of  stations  in  each  0.1  log  amplitude  bin.  Horizontal  axis  is  log 
noise  amplitude  in  Pascals. 

Detection  threshold  maps  were  calculated  for  the  AFTAC  and  LANL  scaling  relations.  The 
results  are  shown  in  Figures  12  to  15.  Figures  12  and  13  were  calculated  with  a  minimum  S/N  of 
2.0,  and  Figures  14  and  15  with  a  minimum  S/N  of  1.5.  For  the  LANL  model  (Figures  12  and 
14),  the  detection  threshold  is  between  0.3  and  2  KT.  For  the  AFTAC  model  (Figures  13  and  15), 
the  detection  threshold  ranges  from  about  0.2  KT  to  about  1  KT.  These  scaling  relations  predict, 
therefore,  that  the  infrasound  detection  threshold  for  the  IMS  network  is  less  than  the  design  goal 
of  one  kiloton  in  most  locations,  but  higher  in  some  regions,  particularly  in  broad  ocean  areas. 

The  network  simulation  results  indicate  that  the  detection  threshold  of  the  future  IMS  infrasound 
network  may  be  somewhat  higher  than  the  one-kiloton  design  goal  in  some  locations.  This  result 
depends,  of  course,  on  a  number  of  assumptions  that  went  into  the  simulations.  We  made  the 
assumptions  that  a  four  element  array  leads  to  a  factor  of  two  improvement  in  signal/noise  ratio, 
that  the  noise  levels  at  all  stations  are  independent,  and  that  a  signal  will  be  identified  with  a 
signal  to  noise  ratio  of  1.5-2.  In  general,  these  assumptions  are  optimistic,  although  experienced 
analysts  may  be  able  to  detect  a  signal  at  lower  S/N  ratios.  Improvements  could  also  be  made  in 
the  signal  and  noise  modeling.  Two  improvements  in  particular  which  would  make  the 
simulations  more  realistic  are: 

1.  including  stratospheric  winds,  which  would  have  the  effect  of  improving  detection  in 
some  directions  and  degrading  it  in  others. 

2.  including  actual  noise  distributions  at  each  station  instead  of  an  average  with  a 
Gaussian  distribution. 
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Figure  12.  Contours  showing  detection  thresholds  with  a  90%  level  of  confidence  for  detection  at  2  infrasound  stations  with  station  dependent  noise  levels  using 
the  LANL  scaling  relation  (Equation  2).  The  intervals  are  logarithmically  spaced  with  labels  in  kilotons.  Minimum  S/N  is  2. 
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Figure  14.  Contours  showing  detection  thresholds  with  a  90%  level  of  confidence  for  detection  at  2  inffasound  stations  with  station  dependent  noise  levels  using 
the  LANL  scaling  relation  (Equation  2).  The  intervals  are  logarithmically  spaced  with  labels  in  kilotons.  Minimum  S/N  is  1.5. 
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Figure  15.  Contours  showing  detection  thresholds  with  a  90%  level  of  confidence  for  detection  at  2  infrasound  stations  with  station  dependent  noise  levels  using 
the  AFT  AC  scaling  relation  (Equation  2).  The  intervals  are  logarithmically  spaced  with  labels  in  kilotons.  Minimum  S/N  is  1.5. 


Figure  16,  for  example,  shows  the  measured  noise  levels  at  the  Los  Alamos  infrasound  station  as 
a  function  of  time  over  a  seven  month  period  starting  January  1999  as  a  time  series  and  as  a 
frequency  distribution.  Two  features  are  apparent.  First,  the  noise  levels  are  time  dependent,  with 
higher  noise  levels  in  winter  than  in  summer.  Second,  the  noise  distribution  is  skewed,  and  much 
higher  than  average  noise  levels  are  common. 


DLI01  Infrasound  measurements 


Figure  16.  Noise  measurements  at  the  Los  Alamos  infrasound  station  for  7  month  period  starting  January,  1999. 
The  top  figure  shows  the  logarithm  of  the  RMS  noise  amplitude  taken  at  noon  each  day  since  the 
beginning  of  1999.  The  bottom  figure  shows  the  frequency  distribution  of  the  noise  measurements.  The 
distribution  of  noise  is  both  skewed  and  time  dependent. 
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Figure  17.  Threshold  contours  of  yield  in  kilotons  for  90%  detection  at  2  or  more  stations  for  the  Gaussian  noise 
model  (left)  and  for  the  non-Gaussian  noise  model  (right).  Simulations  were  performed  with  XNICE. 


NetSim  does  not  have  the  capability  to  model  non-Gaussian  distributions,  but  they  can  be 
modeled  with  XNICE  (Barker,  1996;  Barker  et  al,  1994)  (XNICE  stands  for  the  X  Window 
version  of  the  Network  Identification  Capability  Estimation  system,  which  can  model  network 
identification  as  well  as  detection  capability).  We  modified  XNICE  to  handle  infrasound  as  well 
as  seismic  data,  and  did  a  test  case  to  assess  the  difference  caused  by  a  realistic  noise  distribution 
such  as  that  shown  in  Figure  16.  We  used  a  constant  mean  noise  level  at  all  stations  of  0.1 
Pascals  and  a  minimum  S/N  of  1.5,  and  performed  two  test  cases:  one  with  a  Gaussian 
distribution  with  a  standard  deviation  of  0.64  (derived  from  the  distribution  in  Figure  16),  and 
one  with  the  noise  distribution  modified  to  have  the  same  shape  as  the  distribution  in  Figure  16  at 
all  stations.  Figure  17  shows  the  results  for  South  America.  The  effect  of  including  the  actual 
noise  distribution  is  to  increase  the  threshold  level  by  about  50%,  which  is  a  significant 
difference.  The  IMS  thresholds  should  be  reevaluated  as  actual  noise  data  becomes  available  for 
the  IMS  infrasound  stations. 


22 


Section  3 

Simulations  of  Pipe  Array  Amplitude  and  Phase  Response  and  S/N  Improvement 

3.1  Introduction. 

Pipes  arranged  in  various  configurations,  leading  to  a  pressure  sensor,  are  used  as  spatial  filters  to 
enhance  atmospheric  infrasound  signals  relative  to  ground  level  pressure  perturbations  advected 
by  the  wind.  We  present  a  fast,  accurate,  computationally  straightforward  means  of  calculating 
pipe  configuration  responses,  used  with  appropriate  noise  models,  to  facilitate  the  design  of  pipe 
array  configurations  for  optimal  signal-to-noise  (S/N)  improvement. 

The  numerical  method  is  based  on  acoustic  propagation  in  a  tube,  with  computations  utilizing 
propagator  matrices.  Experiments  verify  our  assumption  regarding  the  physical  basis  for  the 
modeling  and  provide  empirical  measures  of  dispersion  and  attenuation.  We  estimate  the 
amplitude  and  phase  responses  of  two  different  types  of  pipes:  open-ended  pipes  and  closed- 
ended  pipes  with  acoustic  inlets  along  their  lengths,  and  consider  their  relative  merits.  We  then 
estimate  response  curves  for  a  variety  of  pipe  configurations,  made  of  both  pipe  types. 

Temporal  and  spatial  distributions  of  pressure  at  the  ground  are  complicated  and  difficult  to 
predict,  but  are  the  major  factor  in  signal-to-noise  improvement  for  any  pipe  configuration.  We 
use  a  self-similar  noise  distribution,  as  pressure  variations  at  the  frequencies  of  interest  primarily 
result  from  a  cascade  of  energy  from  larger  to  smaller  eddies.  Taylor's  frozen  turbulence 
hypothesis  allows  us  to  produce  noise  time  series  for  each  pipe  inlet  from  the  spatial  noise 
model.  We  present  the  theoretically  achievable  signal-to-noise  improvements,  under  different 
noise  conditions,  for  pipe  configurations  spanning  a  range  of  shapes  and  sizes. 

This  paper  is  arranged  into  four  distinct  segments.  In  the  first,  we  describe  experiments  designed 
to  test  assumptions  regarding  the  physical  basis  for  the  modeling  and  to  provide  empirical 
measures  of  dispersion  and  attenuation.  In  the  second,  we  describe  the  propagator  matrix  method 
we  use  for  numerical  simulations  of  any  configuration.  We  next  discuss  the  noise  model  used. 
Temporal  and  spatial  distributions  of  pressure  at  the  ground  are  complicated  and  difficult  to 
predict.  We  consider  the  impact  on  S/N  improvement  for  pipe  configurations,  of  assumptions 
made  regarding  the  noise  model  and  of  deviations  from  those  assumptions.  Finally,  we  present 
the  results  of  simulations  of  amplitude  and  phase  response  curves  for  two  types  of  pipes  and  of 
arrays  of  such  pipes,  and  of  S/N  improvement  for  the  arrays  under  a  range  of  noise  conditions, 
from  which  we  draw  conclusions  regarding  array  design  and  areas  for  future  work. 

3.2  Background. 

3.2.1  The  Basic  Problem . 

Infrasonic  signals  of  interest  for  monitoring  nuclear  explosions  and  for  many  natural  phenomena 
(e.g.  volcanic  eruptions,  bolides,  aurora,  sonic  booms,  and  microbaroms;  e.g.  Wilson  et  al,  1996) 
generally  have  long  wavelengths,  75  meters  to  35  km,  propagate  at  acoustic  velocities  and  arrive 
at  steep  incidence,  but  overlap  in  frequency  with  pressure  fluctuations  due  to  turbulent  eddies 
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near  the  ground  which  travel  at  the  mean  wind  speed  (e.g.  Daniels,  1959).  Arrays  of  pipes  with 
acoustic  inlets  that  span  a  large  area  and  connect  to  a  central  manifold  where  the  pressure  is 
measured  are  used  to  filter  out  the  spatially  incoherent  pressure  fluctuations  near  the  ground  and 
so  enhance  the  amplitude  of  steeply  incident,  long  wavelength  signals  relative  to  the  noise  (e.g. 
Daniels,  1959;  Burridge,  1971;  Grover,  1971).  This  technique  is  effective,  but  not  yet  thoroughly 
quantified. 

3.2.2.  Previous  Work:  Early  Pipe  Arrays. 

Much  of  the  research  on  pipe  arrays  for  infrasound  instruments  was  performed  prior  to  the  1963 
Limited  Test  Ban  Treaty,  which  banned  nuclear  explosions  in  the  atmosphere,  outer  space,  and 
under  water.  Daniels  (1959)  described  the  use  of  pipe  arrays  to  improve  S/N,  where  the  noise  is 
due  to  short  wavelength  atmospheric  pressure  fluctuations  convected  by  wind.  He  used  a  tapered 
linear  pipe  with  inlets  along  its  length;  their  acoustic  resistance  matched  to  that  of  the  local  pipe 
dimensions  to  prevent  extraneous  reflections  of  signal  within  the  pipe.  Burridge  (1971)  and 
Grover  (1971)  addressed  the  design  and  performance  of  such  arrays.  Burridge  presented  a 
mathematical  model  of  the  response  of  a  linear  pipe,  and  demonstrated  that  non-tapered  pipe 
would  perform  as  well  as  tapered.  He  also  considered  the  performance  of  a  circular  pipe  with 
inlets.  Grover  (1971)  reported  on  the  actual  performance  of  linear  and  circular  pipes  with 
acoustic  inlets.  Both  noted  that  instruments  should  span  sufficient  area  that  noise  coherence  is 
minimized,  but  not  so  much  that  the  phase  is  incoherent  between  signals  received  at  different 
inlets.  Other  design  considerations  they  present  include  the  importance  of  choosing  inlets 
acoustic  resistances  so  that  they  balance  the  contribution  of  nearer  and  more  distant  inlets,  and 
choosing  the  acoustic  impedance  of  the  pipe  for  a  desired  frequency  response. 

3.3  The  Physical  Basis  for  Modeling  a  Pipe  Array  Response. 

Our  goal  is  to  build  on  the  work  described  above  to  determine  the  response  of  pipe  arrays  and  to 
develop  design  criteria  for  optimizing  the  S/N  improvement.  The  first  step  is  to  determine  the 
physical  basis  for  modeling  signal  propagation  within  a  pipe.  The  physical  model  used  by 
Burridge  (1971)  was  that  of  Benade  (1968),  for  acoustic  wave  propagation  in  a  cylindrical 
conduit.  An  alternate  model  however  exists.  For  low  Reynold’s  number  incompressible  fluid 
flow  in  a  cylindrical  conduit,  resistance  to  flow  is  by  viscous  drag  with  zero  velocity  along  the 
boundary.  This  causes  laminar  flow  with  a  parabolic  velocity  gradient  increasing  radially  inward, 
and  is  called  Poiseuille  flow.  Grover  (1971)  demonstrated  that  the  acoustic  resistances  of  inlet 
ports  were  well  approximated  by  Poiseuille  flow,  and  Poiseuille  flow  was  assumed  to  describe 
the  attenuation  in  simulations  of  the  response  and  noise  reducing  capabilities  of  permeable  hose 
by  McLaughlin,  et  al.,  1997. 

We  have  performed  several  experiments  to  determine  how  the  infrasound  signal  propagates 
within  hoses  or  pipes.  The  data  collected  also  allow  us  to  empirically  determine  parameter  values 
for  use  in  simulations. 
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3.3.1  Velocity  and  Dispersion  are  Consistent  with  Acoustic  Wave  Propagation  in  a  Cylinder. 

To  assess  the  validity  and  accuracy  of  both  the  pressure  gradient  driven  flow  and  acoustic  wave 
models,  we  measure  the  velocity  using  two  10-meter  lengths  of  1.6  cm  diameter  permeable  hose 
cut  from  the  same  50  meter  length  (Figure  18).  The  hoses  were  parallel,  within  a  couple 
centimeters  of  each  other,  and  each  was  capped  at  one  end  and  connected  at  the  other  end  to 
identical,  calibrated  pressure  transducers.  A  30-meter  section  of  1.6  cm  diameter  impermeable 
hose  was  placed  between  one  of  the  hoses  and  the  transducer,  allowing  us  to  measure  the  velocity 
and  dispersion  due  to  propagation  through  impermeable  hose  independent  of  any  effects  of 
permeability. 


/ . /I 


30  meters  coiled 
impermeable  hose 


Figure  18.  Layout  of  experiments  described  in  text. 


We  induced  pressure  changes  within  a  narrow,  sealed,  1 1  meter  long  room  and  measured  them 
with  the  apparatus  shown  in  Figure  18,  to  provide  input  signals  between  0.1  and  10  Hz,  with 
amplitudes  from  0.05  to  1  Pascal.  Thus,  these  tests  provide  a  measure  of  the  hose  response  to 
signals  in  the  frequency  and  amplitude  range  of  interest  for  infrasound  monitoring. 

Signals  from  the  two  10  meter  permeable  hoses,  without  the  extra  length  of  impermeable  hose 
shown  in  Figure  18,  correlate  very  well  with  zero  time  lag  between  them  (Figure  19).  When  the 
extra  30  meters  of  impermeable  hose  is  included,  the  signals  also  correlate  well,  but  with  a 
significant  frequency-dependent  time  lag  (Figure  19).  The  velocities  indicate  that  the  infrasound 
signals  propagate  as  acoustic  waves  in  the  hoses,  not  as  pressure  gradient  driven  flow. 
Unfortunately,  significant  deviations  in  signal  amplitude  between  supposedly  identical  10  m 
permeable  hoses  (for  the  experimental  setup  without  the  30  m  impermeable  hose)  preclude 
measurement  of  attenuation  as  a  function  of  frequency.  We  experimented  with  a  few  different 
segments  of  hose  cut  from  a  single  originally  50  m  long  hose.  The  differences  are  presumably 
due  to  variations  in  leak  rate  between  the  hoses.  Such  variability  is  important  to  document  as  it 
necessitates  the  use  of  other  materials  when  the  instrument  response  must  be  precisely  known. 

The  dispersion  is  greater  than  predicted,  with  group  velocities  measured  by  cross-correlation 
approaching  150  m/s  near  0.1  Hz.  Predictions  are  based  on  the  linear  transmission  line  formalism 
used  to  describe  the  acoustical  properties  of  a  tube  (Benade,  1968;  Haak  and  de  Wilde,  1998). 
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The  greater  than  predicted  dispersion  may  have  the  same  cause  as  the  greater  than  predicted 
attenuation,  discussed  next,  as  the  theoretically  dispersion  and  attenuation  are  related  through  the 
propagation  constant  of  the  hose. 


Identical  Hoses  Extra  30m  Hose 


Figure  19.  Segments  of  signals  recorded  by  the  apparatus  shown  in  Figure  18,  filtered  from  1-5  Hz  (top)  and  from 
0.05  -  0.2  Hz  (bottom).  Those  in  the  left  column  were  recorded  with  no  extra  hose,  and  the  maximum 
cross-correlation  of  each  is  at  zero  lag.  Those  on  the  right  had  a  30  m  impermeable  hose  interspersed 
between  the  end  of  one  permeable  hose  and  the  transducer  (dotted  line).  In  that  case,  the  maximum 
cross-correlations  are  at  0.080  (top)  and  0.122  seconds  lag  (bottom),  indicating  significant  dispersion  of 
the  signal. 


3.3.2  Measured  Attenuation  is  much  Greater  than  Predicted. 

The  loss  mechanisms  for  Poiseuille  flow  and  acoustic  waves  in  a  cylinder  are  similar,  in  that  they 
both  lose  energy  in  the  viscous  boundary  layer  along  the  cylinder  wall.  Further,  for  frequencies  of 
interest,  infrasound  wavelengths  are  greater  than  the  longest  pipes  being  considered.  Thus,  all  the 
air  in  the  hose  may  effectively  be  moved  as  a  unit  due  to  the  pressure  pulse  applied  at  one  end 
and  so  attenuation  would  result  from  drag  within  the  boundary  layer  between  the  wall  and  the 
mass  of  air,  similar  to  the  attenuation  for  pressure  gradient  driven  flow.  Finally,  we  note  that  the 
formulas  for  acoustic  propagation  based  on  transmission  line  theory  are  for  an  infinitely  long 
tube.  As  we  discuss  below,  a  steady  state  boundary  layer  may  not  develop  for  acoustic  waves  in 
the  length  and  radii  tubes  considered  here.  Thus,  although  the  velocity  and  dispersion 
measurements  demonstrate  that  the  signal  propagates  as  an  acoustic  wave,  if  the  boundary  layer 
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is  laminar,  the  attenuation  predicted  for  Poiseuille  flow  could  be  more  appropriate  for  these  long 
wavelength  signals. 

We  measured  attenuation  directly  by  recording  pressure  signals  before  and  after  propagation 
along  30,  40,  and  60  m  lengths  of  impermeable  hose,  using  the  setup  illustrated  in  Figure  21.  As 
with  the  velocity  and  dispersion  measurements,  these  results  provide  insight  into  the  controlling 
physical  processes  as  well  as  permitting  empirical  determination  of  parameter  values  for  use  in 
modeling. 


Figure  20.  Theoretical  group  velocity  curves  for  acoustic  waves  in  cylindrical  conduits  (Benade,  1968;  Haak  and  de 
Wilde,  1998),  and  group  velocity  measurements  from  experiments  such  as  illustrated  in  Figure  18.  The 
solid  and  dashed  curves  are,  respectively,  group  velocities  for  acoustic  waves  in  a  cylinder  for  the 
viscosity  of  air  and  for  a  7  times  higher  viscosity  term,  to  account  for  possible  non-laminar  resistance  to 
shear  and  inertial  terms.  The  two  symbols  represent  dispersion  measurements  from  two  distinct 
repetitions  of  the  experiments. 


Figure  21.  Layout  of  experiments  to  measure  attenuation.  We  used  this  setup  with  30,  40,  and  60  meter  lengths  of 
impermeable  hose.  The  extra  30  m  length  of  hose  at  the  end  ensured  that  reflections  from  the  end  did  not 
affect  the  results. 
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The  signal  frequency  typically  peaked  around  2-4  Hz  in  these  experiments.  Figure  22  shows  the 
observed  values  and  the  predicted  attenuation  for  3  Hz  acoustic  waves,  using  two  values  for 
viscosity,  for  various  length  hoses. 


Figure  22.  Predicted  and  observed  signal  decay  as  a  function  of  length  of  the  hoses.  The  solid  line  indicates  the 

decay  predicted  for  acoustic  waves  in  1 .6  cm  diameter  hoses  for  3  Hz  waves.  The  dashed  line  indicates 
the  decay  predicted  using  a  larger  viscosity  value  empirically  determined  to  better  fit  the  data.  The  error 
bars  indicate  2  standard  deviation  uncertainty  bounds  for  observed  attenuation  for  30, 40,  and  60  meter 
long  hoses. 

Neither  attenuation  nor  dispersion  is  modeled  accurately  by  the  equations  for  acoustic 
propagation  in  a  cylinder.  Therefore,  for  numerical  simulations  we  use  an  empirically  determined 
measure  of  viscosity  that  improves  the  data  fit  (Figures  20  and  22). 

These  results  are  consistent  with  those  of  Stevens  et  al.  (1998),  in  which  the  decay  of  Helmholtz 
oscillations  within  tubes  was  much  greater  than  predicted  for  viscous  shear  in  a  laminar  boundary 
layer.  They  also  found  that  the  attenuation  is  independent  of  whether  a  hose  is  permeable  or 
impermeable.  That  is,  the  amplitude  loss  due  to  leakage  through  inlets,  or  permeable  hose,  was 
separable  from  that  due  to  propagation  losses,  and  the  propagation  losses  were  independent  of  the 
conduit  permeability.  The  importance  of  this  is  that  air  movement  through  inlets  does  not 
interfere  with  propagation  through  the  pipes,  so  the  same  attenuation  constant  can  be  used  in 
simulations  of  either  type  of  hoses,  or  for  pipes  with  or  without  inlets. 
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3.3.3  Functional  Dependence  of  Dispersion  and  Attenuation  on  the  Propagation  Constant. 


Attenuation  is  proportional  to  ere(n,  while  phase  velocity  is  given  by  Vp=aVIm(r),  where  T  is  the 
propagation  constant  discussed  in  the  following  section.  While  we  use  the  exact  expressions  for 
T  in  our  calculations,  the  much  simpler  long  period,  small  tube  approximation  for  T  is 
appropriate  for  much  of  the  frequency  band  of  interest  and  is  more  illuminating  for  this 
discussion.  In  that  approximation,  T  approaches 


r-> 


(17) 


where  co  is  angular  frequency,  c  is  sound  velocity  in  air,  y  is  the  ratio  of  specific  heats  at  constant 
pressure  and  volume  Cp/Cv,  and  rv  is  the  ratio  of  the  cylinder  radius  to  the  viscous  boundary  layer 
thickness  (Benade,  1968;  Haak  and  de  Wilde,  1998).  That  ratio  is  given  by 

(18) 
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where  p  and  r)  are  the  density  and  viscosity  of  air,  and  a  is  the  cylinder  radius  (Benade,  1968). 


While  useful  models,  even  the  complicated  exact  expression  for  T  relies  on  assumptions  that  are 
violated  in  practice.  We  discuss  two.  First,  the  expressions  are  for  infinitely  long  tubes,  and  so  do 
not  take  into  account  that  steady  state  may  not  be  reached  in  the  short  pipes  considered  here. 
Thus,  the  equations  lack  inertial  terms  describing  the  acceleration  of  air  in  the  pipes  by  an 
acoustic  wave.  Second,  roughness  in  the  pipe  interior  is  ignored.  Such  roughness  could 
effectively  increase  the  thickness  of  the  boundary  layer.  We  cannot  at  this  time  resolve  the  true 
physics  behind  the  observed  dispersion  and  attenuation.  These  considerations  however  provide 
reason  for  replacing  the  viscosity  term  with  a  larger,  empirically  determined  value.  The  empirical 
term  could  represent  the  additional  resistance  due  to  inertia,  greater  boundary  layer  thickness,  or 
flow  that  is  not  laminar  either  because  steady  state  is  not  reached  or  because  the  wall  roughness 
causes  turbulence.  We  note  that  the  eddy  viscosity,  used  to  describe  viscous-like  effects  in 
turbulent  flow,  can  easily  be  1  to  2  orders  of  magnitude  greater  than  the  fluid  viscosity  (e.g. 
Schlichting,  1960).  Since  T  increases  with  viscosity,  an  increase  in  the  effective  viscosity 
increases  both  the  attenuation  and  dispersion.  This  is  illustrated  in  Figures  20  and  22,  where 
dashed  lines  show  that  both  the  predicted  group  velocity  and  amplitude  loss  are  better  fit  when 
we  incorporate  an  empirical  viscosity  term  7  times  greater  than  that  of  the  viscosity  of  air. 

3.3.4  Summary  and  Conclusions  Regarding  Physical  Model. 

Acoustic  wave  propagation  in  a  cylinder  is  consistent  with  experimental  measurements  of 
velocity.  Attenuation  and  dispersion  however  are  much  greater  than  predicted.  We  therefore  base 
our  simulations  on  the  equations  for  acoustic  propagation,  with  the  viscosity  term  replaced  by  an 
empirically  determined  value  that  provides  numerical  results  that  better  match  observations. 

Different  materials  used  to  construct  spatial  filters  might  produce  quite  different  attenuation,  due 
for  example  to  differences  in  interior  surface  roughness.  Thus,  we  recommend  that  attenuation  be 
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determined  empirically  for  whatever  materials  are  used.  The  theoretical  attenuation  for  Poiseuille 
flow  is  a  likely  minimum  value  for  the  real  attenuation. 

3.4  Numerical  Simulation  Method. 


3.4.1  The  Propagator  Method  for  Sound  Transmission  in  a  Pipe. 


As  developed  and  detailed  by  Burridge  (1971),  a  propagator  method  can  be  used  to  model 
acoustic  propagation  within  a  pipe,  due  to  pressure  fluctuations  at  inlets.  For  a  harmonic  pressure 
fluctuation,  p(x)etca  and  a  harmonic  volume  flux  f(x)e,a*  at  a  distance  x  from  the  left  end  of  a 

(A 

pipe  (Figure  23),  the  resulting  flux-pressure  vector 
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(Burridge,  1971).  T  and  Z0  are  the  propagation  constant  and  characteristic  impedance  of  the  pipe. 

We  generally  use  here  the  expressions  of  Haak  and  de  Wilde  (1998),  which  are  equivalent 
formulations  to  those  of  Benade  (1968).  In  that  formulation,  r  -fzY ,  and  Z0=-Jz/Y ,  where 
Z  =  i(cop/ m2)(l  -  Fve‘^v  )-1  and  Y  =  i{com2  / pc2)(l+(r-l)Fte^') ,  a)  is  angular  frequency,  p  is  density, 
a  is  pipe  radius,  c  is  the  acoustic  velocity  in  air,  and  y  is  the  ratio  of  specific  heats  at  constant 


pressure  and  volume.  Fve = 


2J\(rvfh) 

rvJ-ijo(rvJ-i) 


and  F,e^‘  = 


2Jx(r,^i) 


where  J  are  complex 


Bessel  functions,  rv  is  as  defined  in  equation  2,  and  r,  =  a(copC p  / k)u1  is  the  ratio  of  the  pipe  radius 
to  the  thermal  boundary  layer  thickness,  where  Cp  is  the  specific  heat  of  air  at  constant  pressure, 
and  k  is  the  thermal  conductivity. 


For  an  inlet  at  x,  the  resulting  flux-pressure  vector  propagates  from  one  side  of  the  inlet  (x.)  to 
the  other  side  of  x  (x+)  according  to 
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where  Pe  is  the  pressure  disturbance  at  the  inlet  and  Z*  is  the  impedance  of  the  inlet.  Therefore, 
the  flux-pressure  vectors  at  two  points  along  a  pipe,  xi  and  X2  (Figure  23),  are  related  by  the 
propagator  matrices  representing  the  inlets  and  distances  between  inlets  as  follows, 
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where  Nb  is  the  number  of  inlets  between  x\  and  jc2.  The  inlets  are  numbered  as  1,2,  ...,Nb 
between  x\  and  x2.  The  N’s  are  the  propagator  matrices  at  the  inlets,  defined  in  equation  (4),  M’s 
are  propagator  matrices  between  inlets,  and  1  are  the  propagator  matrices  between  jq  and 

the  first  inlet  and  between  x2  and  the  last  inlet.  A  and  B  are  therefore  the  complete  propagator 
matrices  from  xj  to  x2.  We  next  apply  this  method  to  derive  the  response  of  the  pipe  arrays. 


Figure  23.  Pipe  configuration. 


3.4.2  Formulas  for  Pressure  at  the  Center  of  Close-Ended  Pipe  Arrays. 


For  a  set  of  pipes  arrayed  radially  and  connected  to  a  summing  manifold  at  the  center  (Figure  24, 
upper  left  and  center),  using  the  above  method,  the  pressure  response  at  the  center  is 


P  = 
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(22) 


where  R  is  the  total  number  of  radial  pipes  and  the  A’s  and  B’s  are  the  propagator  matrices  from 
the  pipe  ends  to  the  center  for  each  pipe,  as  defined  in  (19). 

For  a  hexagonal  system  as  shown  in  Figure  24  (upper  right),  the  pressure  response  at  the  center  is 
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where  A(l>  is  the  propagator  matrix  from  the  pipe  end  to  the  center  for  the  i-  radial  impermeable 
pipe,  (l)A(l>,  and  <r>A<l),(r)B<l  are,  respectively,  the  propagator  matrices  from  the  left  and  right 
ends  to  the  center  of  each  transverse  pipe  joined  to  the  radial  pipe  i.  The  right  subscripts  denote 
components  of  matrices. 

Equations  (22)  and  (23)  provide  the  responses  of  two  different  pipe  arrays.  These  are 
straightforward  to  construct  for  any  number  of  radial  and  transverse  pipes,  and  for  irregular 
combinations  of  pipes  with  unevenly  distributed  angles  and  different  radii.  The  responses  at  the 
center  are  represented  as  the  sum  of  all  the  pipe  contributions.  The  method  is  directly  applicable 
as  well  to  configurations  of  open-ended  pipes,  such  as  those  in  the  lower  row  of  Figure  24. 


•  • 


40  m 
120  inlets 


120  inlets 


70m 

144  ports 


Figure  24.  Designs  for  which  we  simulate  S/N  enhancement.  Those  along  the  top  row  are  constructed  of  2.5  cm 
radius  pipes  with  100  acoustic  ohm  inlets  every  meter  (dotted  lines)  connected  to  impermeable  pipe 
(solid  lines).  Pressure  is  measured  at  a  summing  manifold  at  the  center  of  each  configuration.  Designs 
along  the  bottom  row  consist  of  sets  of  essentially  open-ended  impermeable  pipes  arranged  radially  (in 
practice  a  screen  mesh  may  cover  the  ends,  but  this  is  acoustically  insignificant).  The  long  radial  arms 
are  0.95  cm  radius,  and  the  shorter  arms  (and  all  pipes  in  the  18  m  diameter  configuration)  are  0.64  cm 
radius.  Summing  manifolds  are  at  the  center  of  each  set  of  radial  arms,  which  are  then  connected  by  the 
long  radial  pipes  to  a  central  summing  manifold  where  the  pressure  is  measured. 


3.5  Noise  Models. 


3.5.1  Noise  Model  Used. 

The  performance  of  a  pipe  array  depends  on  the  character  of  the  noise,  and  so  meaningful 
modeling  of  infrasound  pipe  array  performance  hinges  on  accurate  knowledge  of  the  temporal 
and  spatial  distribution  of  pressure  fluctuations  over  the  array.  We  discuss  the  physical  basis  for 
the  noise  model  chosen  and  limitations  imposed  by  necessary  simplifications. 

Pressure  variations  across  the  ground  are  assumed  to  be  due  to  and  linearly  related  to  turbulent 
eddies,  which  are  generally  represented  as  perturbations  to  the  mean  wind  velocity  (e.g.  Stull, 
1988;  Kaimal  and  Finnigan,  1994).  The  scale  of  eddies  depends  on  the  boundary  layer  thickness, 
and  within  a  canopy,  depends  on  the  canopy  height.  Large  turbulent  eddies  are  generated  by 
buoyancy  driven  convection,  typically  due  to  heating  of  the  surface  layer  over  land.  Smaller 
eddies  are  then  generated,  as  wind  shear  in  the  boundary  layer  breaks  up  the  large  eddies.  Eddies 
in  the  scale  length  of  interest  for  infrasound  are  in  the  inertial  subrange  of  the  turbulent  kinetic 
energy  (TKE)  spectrum.  That  is,  they  are  due  to  a  cascade  of  energy  from  larger  to  smaller 
eddies.  We  use  a  zeroth  order  Von  Karmann,  or  self-similar,  distribution  to  represent  that 
distribution.  That  is  defined  by  a  correlation  function  of  Ko(r/a),  where  Ko  is  the  zeroth  order 
Bessel  function  (e.g.  Frankel  and  Clayton,  1986). 

Taylor's  frozen  turbulence  hypothesis  is  almost  universally  assumed  in  studies  of  wind  near  the 
ground  surface  (e.g.  Kaimal  and  Finnigan,  1994;  Stull,  1988).  It  assumes  that  eddies  evolve 
slowly  compared  with  the  time  it  takes  for  the  mean  wind  velocity  to  transport  them  past  a  point 
on  the  ground  and  that  all  sizes  of  eddies  move  at  the  same  (wind)  velocity.  The  assumption  is 
ubiquitous  because  it  permits  the  transformation  of  a  time  series  measured  at  one  point  to  an 
estimate  of  the  spatial  dimensions  of  eddies.  Taylor's  hypothesis  holds  to  first  order  in  many 
situations,  and  we  assume  it  for  our  simulations,  which  allows  us  to  transform  the  simulated, 
self-similar,  spatial  pressure  variations  into  time  series  at  each  instrument  port. 

3.5.2  Complexities  in  Real  Noise:  Impact  on  Results  and  Areas  for  Further  Research. 

A  Von  Karmann  noise  model  represents  a  simplification  of  the  real  physics.  Turbulence  in  the 
atmospheric  boundary  layer,  and  particularly  near  the  ground  and  in  canopies  is  very  complex. 
We  consider  briefly  the  potentially  most  important  deviations  from  the  simple  model  we  assume, 
since  complexities  not  incorporated  into  the  noise  model  may  effect  the  accuracy  of  simulations 
of  array  performance. 

Stull  (1988)  reports  that  Taylor's  hypothesis  generally  holds  when  the  turbulent  intensity,  given 
by  the  standard  deviation  of  wind  speed,  is  much  less  than  the  mean  wind  speed.  In  practice,  this 
means  it  breaks  down  under  very  gusty  conditions.  Evolution  of  eddies  over  time  periods  less 
than  it  takes  for  them  to  pass  over  a  point  should  lead  to  less  spatial  coherence,  and  so  better 
performance.  Since  such  conditions  tend  to  exist  at  higher  wind  velocities,  we  might  expect 
performance  to  deteriorate  less  than  linearly  with  the  amplitude  increase  that  comes  with  a  mean 
wind  speed  increase  (e.g.  McDonald,  et  al.,  1971). 
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Kaimal  and  Finnigan  (1994)  note  that  larger  eddies  tend  to  travel  more  rapidly  than  smaller  ones, 
especially  in  plant  canopies.  Taylor’s  hypothesis  assumes  that  all  eddies  move  at  the  same  speed. 
The  impact  of  eddy  size  and  wind  speed  is  not  entirely  intuitive.  It  is  obvious  that  a  large  eddy 
covering  a  significant  portion  of  the  array  will  contribute  to  shorter  period  noise  the  more  rapidly 
it  moves.  Small  eddies  however,  can  also  lead  to  surprisingly  significant  noise  coherence,  and  so 
little  S/N  improvement,  as  the  same  signal  traverses  subsets  of  inlets  at  some  time  lag.  The  size 
and  velocity  together,  of  the  eddies,  will  determine  which  frequencies  are  affected.  If  accurate 
measures  of  such  features  of  the  noise  distribution  are  available  for  a  site,  they  could  be 
incorporated  into  simulations  to  indicate  how  the  S/N  will  be  affected. 

The  TKE  spectrum  may  be  spiked  at  wavenumbers  corresponding  to  particular  obstructions  to 
the  wind.  In  a  canopy,  tree  trunks,  brush,  or  grasses  may  break  down  larger  eddies  into  eddies  too 
small  to  be  at  all  coherent  between  inlets.  Thus,  denser  canopies  with  more  small-scale  structures 
near  the  ground  should  provide  advantages  over  more  open  canopies.  In  the  case  of  larger, 
topographic  obstructions,  Kaimal  (1998)  reports  that  a  “hill  wake  may  be  detectable  many  tens  of 
hill  lengths  downwind.”  This  is  an  important  consideration  in  siting. 

Another  feature  of  wind  in  canopies  that  may  be  important  for  pipe  array  design  is  the  existence 
of  horizontal  and  vertical  gusts  on  the  scale  of  the  canopy  height.  This  produces  a  non-trivial 
third  moment,  or  skewness,  in  the  TKE  spectrum  (e.g.  Kaimal  and  Finnigan,  1994;  Finnigan  and 
Brunet,  1995).  A  vertical  gust  that  spans  a  single  pipe  configuration  could  be  indistinguishable 
from  an  infrasonic  signal  at  that  single  receiver. 

Finally,  we  note  that  the  scale  of  spatial  coherence  of  wind  can  be  many  times  greater  parallel  to 
the  mean  wind  velocity,  than  normal  to  it  (McDonald  et  al,  1971;  Finnigan  and  Brunet,  1995). 

For  sites  Where  the  wind  comes  from  a  single  predominant  direction  for  the  vast  majority  of  the 
time,  better  performance  may  be  achievable  with  fewer  inlets  using  an  asymmetric  array  that 
takes  advantage  of  the  asymmetry  of  the  wind  coherence.  Where  instrument  siting  options  are 
limited,  simulations  that  include  such  complications  in  the  noise  model  will  provide  more 
reliable  predictions  of  performance. 

Because  S/N  improvement  depends  on  decorrelation  over  time  at  a  single  inlet,  and  spatially  over 
multiple  inlets,  the  unmodeled  complexities  of  the  noise  described  above  will  affect  the  accuracy 
of  the  S/N  improvement  simulation.  Thus,  the  discussion  above  should  provide  direction  for 
future  work  with  the  goal  of  enabling  better  modeling  of  the  noise  structure,  including  better 
understanding  its  genesis.  The  first  step  in  that  effort  should  be  to  better  quantify  the  spatial 
statistics  of  pressure  variations,  to  provide  real  spatial  and  temporal  noise  models  with 
correlations  accurately  determined  for  different  types  of  sites.  In  addition,  current  site  surveys  are 
performed  with  a  single  sensor,  requiring  Taylor's  hypothesis  for  transformation  to  produce  a 
spatial  noise  model.  Consideration  should  be  given  to  performing  future  site  surveys  with  an 
array  of  sensors,  in  order  to  capture  rather  than  assume  the  spatial  distribution  of  noise. 
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3.6  Simulations. 


3.6.1  Overview. 

We  perform  simulations  for  the  6  different  configurations  of  pipe  arrays  shown  in  24.  Each  of  the 
basic  configurations  have  been  deployed  by  infrasound  researchers  or  considered  for  deployment 
as  part  of  the  International  Monitoring  System  (Christie,  1999).  These  configurations  are  of  two 
distinct  types.  One  type  is  made  up  of  impermeable  pipes  and  close-ended  pipes  with  small  inlets 
of  high  acoustic  impedance  along  their  lengths.  These  are  similar  to  arrays  constructed  of 
permeable  hose.  The  other  type  uses  impermeable  pipes  with  open  ends,  providing  almost  no 
acoustic  impedance  mismatch  with  the  atmosphere.  Both  types  of  configurations  use  one  or  more 
summing  manifolds  where  multiple  pipes  meet.  We  assume  that  the  manifold  volumes  are 
minimal.  If  that  straightforward  engineering  constraint  is  not  met,  reflections  off  the  manifolds 
could  lead  to  poorer  performance. 

3.6.2  Amplitude  and  Phase  Responses  of  Open  and  Closed  End  Pipes. 

The  responses  of  open-ended  pipes  vs.  closed-ended  pipes  with  acoustic  inlets  differ  in  some 
important  respects  (Figures  25  to  29  We  illustrate  these  differences  by  comparing  the  responses 
of  two  pipes,  both  60  meters  long.  One  is  open-ended,  but  with  no  acoustic  inlets.  The  second  is 
closed-ended,  with  one  acoustic  inlet  of  100  acoustic  ohms  resistance  in  the  center.  The 
amplitude  response  of  the  open  pipe  is  flat  up  to  the  frequency  where  resonance  occurs,  while  the 
closed-ended  pipe  acts  like  a  low  pass  filter  (Figure  25  The  comer  of  the  filter  depends  on  the 
acoustic  impedance  of  the  inlets  and  the  radius  of  the  pipe. 

Resonance  interferes  with  high  frequency  signals.  The  first  resonance  peak  for  open-ended  pipes 
occurs  at  a  frequency  of  f=Vaa/2L,  where  L  is  the  pipe  length  and  Va  is  the  acoustic  wave  speed. 
For  closed-ended  pipes  there  is  an  interference  null  at  that  frequency.  Our  highest  frequency  of 
interest  determines  how  long  of  pipes  we  may  use. 


Figure  25.  Amplitude  response  of  a  60  m  open-ended  pipe  (dashed)  and  a  60  m  closed-ended  pipe  with  a  100 
acoustic  ohm  inlet  at  its  center  (solid),  both  with  a  radius  of  1.25  cm. 
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For  a  constant  inlet  impedance,  the  amplitude  response  of  the  closed-ended  pipe  also  varies  with 
radius  (Figure  26).  The  wider  the  pipe  is,  the  smaller  the  amplitude  response  is.  For  a  wide  pipe, 
the  impedance  is  small  compared  to  that  of  the  inlet,  and  so  the  pressure  perturbation  that  can  be 
achieved  within  the  pipe  is  small.  A  physically  intuitive  explanation  for  this  is  that  only  a  small 
pressure  perturbation  is  permitted  by  the  small  inlet,  and  the  larger  the  volume  of  the  pipe,  the 
more  that  signal  is  diluted.  This  is  not  a  problem  for  open-ended  pipes. 


Figure  26.  Amplitude  response  of  60  m  long  closed-ended  pipes  of  four  different  radii,  each  with  a  single  100 

acoustic  ohm  inlet  in  the  pipe  center.  The  first  resonance  peaks  are  slightly  shifted  because  the  acoustic 
velocity  is  less  in  the  narrower  pipe. 

The  phase  response  is  nearly  constant  for  the  widest  close-ended  pipe  because  the  impedance  of 
the  pipe  is  quite  small  for  all  frequencies  of  interest,  relative  to  that  of  the  inlet  (Figure  27).  The 
impedance,  however,  increases  for  narrower  pipes  and  lower  frequencies.  For  sufficiently  narrow 
pipes  and  low  frequency,  the  pipe  impedance  is  similar  to  that  of  the  inlet  and  so  there  is  less 
delay.  Because  there  is  still  a  mismatch  between  the  inlets  and  the  narrow  pipe  response  at  higher 
frequencies,  there  is  a  large  change  in  phase  delay  from  low  to  high  frequency.  Our  calculations 
assume  a  constant  impedance  inlet,  an  approximation  valid  as  long  as  the  impedance  of  the  inlet 
is  much  less  than  that  of  the  pipe.  Thus,  the  approach  of  the  phase  response  curves  to  zero  really 
represents  an  upper  bound  on  the  rate  at  which  they  will  asymptotically  approach  zero  delay.  The 
important  points  demonstrated  here  are  the  trade-offs  between  pipe  radius,  inlet  impedance,  and 
amplitude  and  phase  response.  There  is  also  a  trade-off  between  pipe  radius  and  the  useful 
frequency  band  near  the  first  resonance  peak.  That  is,  at  the  high  frequency  end  of  the  curves,  the 
phase  response  changes  abruptly  at  the  resonance  peak.  The  curve  of  the  narrowest  (lowest 
impedance)  pipe  is  much  smoother.  Thus,  the  same  resonance  peak  affects  a  wider  frequency 
range  in  a  narrower  pipe. 
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100 


Figure  27.  Phase  response  of  four  60  m  long  closed-ended  pipes  of  different  radii,  each  with  a  100  acoustic  ohm 
inlet  at  their  centers. 

The  amplitude  (Figure  28)  and  phase  (Figure  29)  responses  of  the  open-ended  pipes  are  quite  flat 
up  to  the  resonance  peaks. 


Figure  28.  Amplitude  responses  of  60  m  open-ended  pipes  of  various  radii. 
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Figure  29.  Phase  responses  of  60  m  open-ended  pipes  of  various  radii. 


Open-ended  pipes  appear  to  have  the  advantage  over  closed-ended  pipes  with  small  acoustic 
inlets,  as  they  provide  flat  phase  and  amplitude  responses,  although  the  addition  of  multiple  inlets 
does  flatten  the  amplitude  response  of  the  close-ended  designs  somewhat.  Infrasound  spatial 
filters  made  of  pipes  with  inlets  or  permeable  hoses  are  commonly  used  and  have  a  long  history, 
although  the  advantages  of  open-ended  pipes  have  begun  to  be  recognized  (e.g.  Christie,  1999). 
We  simulate  the  performance  of  configurations  with  both  types  of  pipes  in  the  next  section. 

3.6.3  Pipe  Array  Configurations. 

Three  closed-ended  configurations  modeled  include  the  common  radial  arm  design  (e.g.  Noel  and 
Whitaker,  1991),  a  modified  radial  arm  to  prevent  coherent  noise  near  the  center  from 
dominating,  and  a  hexagonal  design  (Christie,  1999)  intended  for  more  widely  distributed  spatial 
sampling  (Figure  24).  The  open-ended  designs  are  variations  on  one  suggested  by  Alcoverro 
(1998),  and  include  a  small  aperture  design  for  low  wind  conditions,  and  two  70  m  diameter 
configurations  with  96  and  144  ports,  intended  for  long  -period  recording  under  high  wind 
conditions  (Christie,  1999). 

3.6.4  Amplitude  and  Phase  Responses  of  Different  Configurations. 

Figure  30  shows  the  amplitude  responses  of  the  configurations  shown  in  Figure  24.  The 
configurations  using  closed-ended  pipes  (hexagonal  and  two  radial  configurations,  light  lines)  all 
act  as  lowpass  filters,  while  the  open-ended  pipe  configurations  have  flat  amplitude  responses  up 
to  the  resonance  frequencies. 
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The  phase  response  of  the  open-ended  pipes  is  flat  up  to  the  resonance  frequency.  The  influence 
of  the  resonance  frequency  however  extends  quite  a  ways.  The  extent  of  that  influence  could  be 
minimized  without  badly  effecting  the  amplitude  response,  by  the  use  of  slightly  wider  radius 
pipe  (see  Figure  27). 


The  closed-ended  pipes  have  phase  delays  of  approximately  50°  at  1  Hz,  well  before  the 
resonance  frequency  has  any  effect.  The  phase  delay  however  is  not  as  severe  as  that  for  a  single 
inlet  (Figure  27),  since  the  impedance  mismatch  between  the  sum  of  the  120  inlets  and  the  pipe  is 
less  than  it  is  for  a  single  inlet.  These  configurations  consist  of  2.5  cm  radius  pipes.  Narrower 
pipes  would  reduce  the  phase  delay  further  and  provide  an  overall  higher  amplitude  response  for 
these  closed-pipe  configurations. 


Frequency 


Figure  31.  Phase  response  of  the  configurations  shown  in  Figure  24. 
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3.6.5  Signal  to  Noise  Improvement. 

We  compare  the  improvement  in  S/N  ratio  measured  at  the  central  manifold  of  each 
configuration  shown  in  Figure  24  over  that  recorded  with  no  pipe  array.  If  the  noise  were 
spatially  uncorrelated,  the  S/N  improvement  would  simply  be  Vn  ,  where  n  is  the  number  of  ports 
or  inlets.  We  use  a  vertically  incident  plane  wave  signal  and  Von  Karmann  distributed  noise  with 
correlation  lengths  of  1  and  10  meters  (Figure  32).  As  discussed  above,  we  assume  that  the  noise 
field  remains  stable  as  it  is  carried  passed  the  instrument.  This  feature  of  the  noise  has  a 
tremendous  impact  on  the  performance  of  any  instrument. 

Both  greater  correlation  lengths  and  lower  wind  speeds  lead  to  lower  frequency  pressure 
variations  at  any  single  spatial  position  (Figure  33).  This  is  due  to  the  advection  of  stable 
correlated  noise  and  will  impact  each  instrument  similarly. 

Advection  of  spatially  correlated  noise  leads  to  temporal  correlations  at  individual  ports  or  inlets. 
That  is,  pressure  variations  at  two  ports  aligned  along  the  wind  direction  will  be  highly  correlated 
at  a  time  lag  of  x  =  D/Vw,  where  D  is  the  distance  between  ports,  and  Vw  is  the  wind  speed.  At 
signal  periods  much  greater  than  t,  the  two  ports  have  highly  correlated  inputs  and  will  not 
contribute  to  increasing  the  S/N.  Ports  separated  by  less  than  the  correlation  length,  in  any 
direction,  will  also  have  highly  correlated  inputs  at  zero  time  lag.  Thus,  we  see  that  smaller 
diameter  configurations,  longer  correlation  length  noise,  and  higher  wind  speeds  lead  to  poorer 
low  frequency  performance  (Figure  34).  In  the  extreme  case  of  noise  with  a  correlation  length 
much  greater  than  the  array  diameter,  the  noise  would  be  coherent  over  much  of  a  pipe  array  and 
would  sum  as  coherently  as  steeply  incident  signals. 

On  the  other  hand,  at  frequencies  of /=  Vw/2-D,  the  signals  will  be  negatively  correlated  at  ports 
aligned  with  the  wind  direction  and  separated  by  a  distance  D,  as  defined  above.  This  occurs  in 
some  of  the  simulations,  where,  at  certain  frequencies  we  observe  a  better  than  Vn  improvement 
in  S/N,  where  n  is  the  number  of  ports  (Figure  34).  Even  small  deviations  from  our  ideal  noise 
model  will  make  that  less  likely  to  occur  in  real  noise  conditions.  The  spectra  shown  are  good 
guides  to  expected  performance,  but  such  fine  details  should  not  be  over-interpreted. 

We  do  not  complicate  the  presentation  of  results  by  scaling  the  noise  amplitude  by  the  wind 
speed,  although  higher  wind  speeds  are  clearly  associated  with  higher  amplitude  pressure 
variations  (e.g.  McDonald  et  al,  1971).  Such  complications  could  be  simply  factored  into  the 
results  after  calculations  for  unit  signal  and  noise  levels  have  been  made,  so  comparison  can  be 
made  of  effective  S/N  ratios  under  different  wind  conditions. 
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Figure  32.  Seventy-five  by  seventy-five  meter  samples  of  spatially  uncorrelated,  and  1  and  10  meter  correlation 
length,  Von  Karmann  distributed  noise. 


10  m  Correlation  Length,  1  m/s  10  m  Correlation  Length,  6  m/s 


Time  (seconds)  Time  (seconds) 


Figure  33.  Pressure  variation  time  series  observed  at  a  single  point  due  to  advection  of  the  three  different  noise 
fields  (Figure  32)  at  1  m/s  (left  column)  and  at  6  m/s  (right  column). 
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Figure  34.  S/N  improvement  for  six  different  configurations  (Figure  24),  for  1,  and  10  m  correlation  length  self¬ 
similar  noise,  carried  past  the  instruments  by  1  and  6  m/s  winds.  Maximum  radius  spanned  and  number 
of  ports  (for  open-ended  pipes)  are  listed.  Each  of  the  close-ended  pipe  configurations  has  120  inlets. 
For  uncorrelated  noise,  S/N  improvement  will  simply  be  a  constant  Vn  improvement,  where  n  is  the 
number  of  ports. 


3.7  Conclusions 

Acoustic  propagation  in  a  cylindrical  conduit  provides  a  good  physical  basis  for  modeling  the 
response  to  an  infrasound  signal  of  an  array  of  pipes.  Actual  attenuation  and  dispersion  however 
can  be  much  greater  than  predicted,  indicating  that  the  particular  pipes  used  should  be  calibrated. 

A  self-similar  noise  distribution  provides  an  appropriate  and  simple  noise  model,  as  pressure 
variations  in  the  frequency  of  interest  primarily  result  from  a  cascade  of  energy  from  larger  to 
smaller  eddies.  Taylor's  frozen  turbulence  hypothesis  allows  us  to  produce  noise  time  series  at 
each  pipe  inlet  from  the  spatial  noise  model.  Both  assumptions  together  enable  us  to  produce 
noise  time  series  at  each  pipe  inlet  that  have  realistic  relationships  with  the  noise  series  at  all  the 
other  inlets. 
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Open-ended  pipes  hold  most  of  the  advantages  over  closed-ended  pipes  with  high  impedance 
acoustic  inlets.  Specifically,  the  open-ended  pipes  have  essentially  flat  amplitude  responses, 
while  the  closed-ended  pipes  act  as  lowpass  filters.  The  open-ended  pipes  also  have  flat  phase 
responses,  while  the  phase  responses  of  the  closed-ended  pipes  vary  much  more  across  the 
spectrum  of  interest.  In  general,  the  acoustic  impedance  of  inlets  and  radii  of  closed-end  pipes 
must  be  carefully  matched  to  achieve  a  desired  response  curve,  and  even  then,  it  will  not  be  flat. 

The  advantages  described  above  suggest  that  for  most  purposes,  open-ended  pipes  will  be 
preferable.  For  permanent  installations,  pipes  are  certainly  preferable,  but  permeable  hose  has 
been  a  mainstay  of  temporary  field  campaigns  for  volcanic  monitoring.  Our  results  indicate  that 
flatter  amplitude  and  phase  response  functions  could  be  obtained,  while  maintaining  the  low 
expense  and  simple  transport  and  set-up  of  hoses,  by  switching  from  permeable  to  open-ended 
impermeable  hoses. 

The  most  important  design  criterion  for  enhancing  S/N  is  having  sufficient  spatial  sampling. 
Extra  ports,  however,  provide  no  advantage  if  spaced  more  closely  than  the  noise  correlation 
length.  All  configurations  performed  well,  and  similarly,  at  high  frequencies,  as  the  high 
frequency  noise  has  the  least  spatial  correlation.  The  differences  in  performance  are  very 
significant  at  the  lower  frequencies.  The  two  70  m  diameter  configurations  provide  the  best 
performance  in  our  simulations  of  S/N  enhancement,  and  their  performances  were  nearly 
identical  even  though  one  had  144  ports  and  the  other  had  only  96.  The  18  m  diameter 
configuration  performed  the  worst  at  low  frequency. 

We  have  developed  a  rapid,  accurate,  method  of  estimating  the  performance  of  any  pipe 
configuration  used  for  infrasound  recording.  The  method  is  based  on  acoustic  propagation  in  a 
cylinder,  the  computations  are  based  on  propagator  matrices,  and  results  match  experimental 
data.  This  tool  will  permit  numerical  simulation-based  optimization  of  new  designs  before 
deployment,  thus  speeding  up  development  and  testing  of  new  designs.  It  will  also  enable 
optimization  of  configurations  where  noise  conditions  are  well  known,  and  estimation  of 
response  curves  for  existing  instruments.  The  most  important  step  that  can  be  taken  now  to 
improve  instrument  performance  is  to  obtain,  and  utilize  in  such  modeling,  more  complete  and 
accurate  noise  models. 
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Section  4 

Electronic  Addendums 

There  are  two  electronic  addendums  to  this  report.  They  are: 

1.  The  complete  digitized  data  set  of  infrasound  recordings;  and 

2.  The  final  report  from  IDG  which  includes  additional  measurements  derived  from  the  data  and 
discussion  of  azimuthal  variations  due  to  wind. 


Section  5 
References 


Alcoverro,  B.  (1998),  Acoustic  filters  design  and  experimental  results,  Proceeding:  Workshop  on 
Infrasound,  Commissarriat  a  l'Energie  Atomique,  Bruyeres  le  Chatel,  France.  (Unclassified) 

Barker,  T.G.  (1996),  “Xnice:  A  System  for  Assessing  Network  Identification  Performance,” 
Maxwell  Laboratories  Scientific  Report  to  Phillips  Laboratory,  PL-TR-96-2087.  (Unclassified) 

Barker,  T.  G.,  K.  L.  McLaughlin,  and  J.  L.  Stevens  (1994),  "Network  Identification  Capability 
Evaluation  (NICE)",  S-CUBED  Technical  Report  to  the  Advanced  Research  Projects  Agency 
SSS-TR-94-14701,  July.  (Unclassified) 

Benade,  A.H.  (1968),  On  the  propagation  of  sound  waves  in  a  cylindrical  conduit,  J.  Acoust.  Soc. 
Am.,  44,  616-623.  (Unclassified) 

Blandford,  R.  R.  and  D.  A.  Clauter  (1995),  “Capability  Estimation  of  Infrasound  Networks”, 
AFTAC  report.  (Unclassified) 

Blandford,  R.  R.,  D.  A.  Clauter,  R.  W.  Whitaker,  and  T.  Armstrong  (1995),  “Infrasound  Network 
Options,”  AFTAC/DOE  White  Paper.  (Unclassified) 

Brown,  D.  J.  (1999)  “Summary  of  Infrasound  Source  Location  Meeting:  San  Diego,  November 
9-10,  1998,  Center  for  Monitoring  Research  Technical  Report  CMR-99/02,  January. 
(Unclassified) 

Burridge,  R.  (1971),  On  the  propagation  of  sound  waves  in  a  cylindrical  conduit,  Geophys.  J.  R. 
Astr.  Soc.,  26,  53-69.  (Unclassified) 

Christie,  D.R.  (1999),  Wind-noise  reducing  pipe  arrays,  memo  of  the  Provisional  Technical 
Secretariat,  CTBTO,  Vienna,  Austria.  (Unclassified) 

Clauter,  D.  A.  and  R.  R.  Blandford  (1998),  “Capability  Modeling  of  the  Proposed  International 
Monitoring  System  60-Station  Infrasonic  Network,”  proceedings  of  the  Infrasound  workshop  for 
CTBT  monitoring,  Santa  Fe,  New  Mexico,  August  25-28,  1997,  LANL  report  number  LA-UR- 
98-56.  (Unclassified) 

Daniels,  F.B.  (1959),  Noise  reducing  line  microphone  for  frequencies  below  1  c/s,  J.  Acoust.  Soc. 
Am.,  31,  529.  (Unclassified) 

Finnigan,  J.J.  and  Y.  Brunet  (1995),  Turbulent  airflow  in  forests  on  flat  and  hilly  terrain,  In  Wind 
and  Trees,  pp  3-40,  M.P.  Coutts  and  J.Grace  (eds.),  Cambridge  University  Press.  (Unclassified) 


45 


Frankel,  A.  and  R.  Clayton  (1986),  Finite  difference  simulations  of  seismic  scattering: 
Implications  for  the  propagation  of  short-period  seismic  waves  in  the  crust  and  models  of  crustal 
heterogeneity,  J.  Geophys.  Res.,  91,  6465-6489.  (Unclassified) 

Grover,  F.H.  (1971),  Experimental  noise  reducers  for  an  active  microbarograph  array,  Geophys. 

J.  R.  Astr.  Soc.,  26, 41-52.  (Unclassified) 

Haak  H.W.  and  G.J.  de  Wilde,  On  the  propagation  of  low  frequency  sound  waves  in  porous 
hoses  used  in  infrasound  detection  systems,  pre-print.  (Unclassified) 

Kaimal,  J.C.  and  J.J.  Finnigan  (1994),  Atmospheric  Boundary  Layer  Flows:  Their  Structure  and 
Measurement,  Oxford  Univ.  Press.  (Unclassified) 

McDonald,  J.A.,  EJ.  Douze,  and  Eugene  Herrin  (1971),  The  structure  of  atmospheric  turbulence 
and  its  application  to  the  design  of  pipe  arrays,  Geophys.  J.  R.  Astr.  Soc.,  26,  99-109. 
(Unclassified) 

McKisic,  J.  M.  (1997),  “Infrasound  and  the  Infrasonic  Monitoring  of  Atmospheric  Nuclear 
Explosions,”  Tracor  Applied  Sciences  final  report  to  Phillips  Laboratory,  PL-TR-97-2123, 
February.  (Unclassified) 

McLaughlin,  K.L.,  B.  Shkoller,  D.E.  Wilkins,  and  T.G.  Barker  (1997),  Advances  in  numerical 
methods.  Maxwell  Technologies  Report  MFD-TR-97-15870.  (Unclassified) 

Pierce,  A.  D.  and  J.  W.  Posey  (1971),  ‘Theory  of  the  Excitation  and  Propagation  of  Lamb’s 
Atmospheric  Edge  Mode  from  Nuclear  Explosions,”  Geophys.  J.  R.  astr.  Soc.,  26,  341-368. 
(Unclassified) 

Posey,  J.  W.  and  A.  D.  Pierce  (1971),  “Estimation  of  nuclear  explosion  energies  from 
microbarograph  records,”  Nature,  232,  253.  (Unclassified) 

Noel,  S.D.  and  Whitaker,  R.W.  (1991),  Comparison  of  noise  reduction  systems,  LANL  report 
LA- 12008-MS  UC-700.  (Unclassified) 

Schlichting,  H.  (1960),  Boundary  layer  theory,  Publ.  McGraw  Hill.  (Unclassified) 

Sereno,  T.  J.,  S.  R.  Bratt,  and  G.  Yee  (1990),  “NetSim:  A  computer  program  for  simulating 
detection  and  location  capability  of  regional  seismic  networks,”  SAIC  Annual  Technical  Report 
to  DARPA,  SAIC  90/1 163,  March.  (Unclassified) 

Stevens,  J.L.,  D.  Adams,  G.E.  Baker,  H.  Xu,  J.  Morris,  I.  Divnov,  V.  Bourchik,  and  I.  Kitov, 
(1999),  Infrasound  scaling  and  sttenuation  relations  from  Soviet  explosion  data  and  instrument 
design  criteria  from  experiments  and  simulations,  Proceedings  21th  Annual  Seismic  Research 
Symposium  on  monitoring  a  CTBT,  185-194.  (Unclassified) 


46 


Stevens,  J.L.,  G.E.  Baker,  J.R.  Murphy,  and  J.  Morris  (1998),  Infrasound  excitation  and 
propagation  research,  Proceedings  2Cfh  Annual  Seismic  Research  Symposium  on  monitoring  a 
CTBT.  (Unclassified) 


Stull,  R  (1988),  An  Introduction  to  Boundary  Layer  Meteorology,  Publ.  Kluwer  Acad.  Press. 
Wexler,  H.  and  W.  A.  Hass  (1962),  “Global  atmospheric  pressure  effects  of  the  October  30, 1961 
explosion,”  J.  Geophys.  Res.,  v.  67,  pp.  3875-3887.  (Unclassified) 

Whitaker,  R.  W.  (1995),  “Infrasonic  Monitoring,”  proceedings  of  the  17th  annual  Seismic 
Research  Symposium  in  Scottsdale,  AZ,  September  12-155,  1995,  997-1000.  (Unclassified) 

Wilson,  R.W.,  J.V.  Olson,  and  R.  Richards  (1996),  Library  of  typical  infrasound  signals,  ENSCO 
Report.  (Unclassified) 


47 


Appendix  A 

Time  constants  of  low  -  (xO  and  high-pass  (Xh)  filters  of  microbarographs. 

Instrument  response  correction.  The  microbarograph  is  a  differential  manomenter  made  up  of 
two  volumes  separated  by  a  diaphragm.  The  first  volume  is  directly  connected  to  atmosphere,  the 
second  one  through  a  high  acoustic  resistance  (capillary  tube).  The  diaphragm’s  flexure  is 
proportional  to  5P.  This  system  is  an  acoustical  high-pass  filter  with  time  constant: 

xh  =  8r)LV  /  7tr4p0, 

where  r)  -  dynamic  viscosity  of  air,  L  and  r  -  length  and  radius  of  capillary  tube,  Po  -  average 
atmospheric  pressure.  High-frequency  (f  »  l/27tXh)  pressure  variations  pass  such  filter  with 
minimal  distortion,  but  low-frequency  pressure  variations  attenuate  and  lead  atmospheric 
pressure  changes  in  phase. 

The  amplitude  Ah(f)  and  phase  (ph(f)  responses  of  the  microbarograph  as  an  acoustical  high-pass 
filter  can  be  written  as: 


_ 1 

V1  +  1  /(27rfxh)T  ’ 


cph  =  arctg  (1  /  2n  f  xh). 


To  restrict  the  frequency  band,  the  microbarographs  (except  MBS)  were  equipped  with  built-in 
electrical  low-pass  filters,  which  attenuate  and  lag  in  phase  the  high  frequency  signal  relative  to 
the  input  signal,  but  allows  to  pass  low  (f  «  1/2txTi)  frequencies.  Amplitude  Ai  and  phase  cpi 
responses  of  the  low-pass  filter  are: 


Ai  =  .  ,  .  ,  <Pi  =  -  arctg  (2n  f  t,). 

^  1  +  (27tfx,) 

The  resulting  responses  are  defined  from  (2)  and  (3)  as: 
A  =  AhX  Ai ,  cp  =  cph  +  cpi ,  A  =  A-ei<p,  j  -  imaginary  unit. 
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Time  constants  ih  and  T|  are  listed  in  the  table  below  for  each  station,  and  are  listed  for  each 
waveform  in  appendix  B. 


Station 

No 

Name 

Device  type 

Tl,S 

Th,S 

Records* 

number1) 

In  all  on  this 
station11 

1 

Dubna 

MBS 

0 

32.2 

2 

5 

LMB-P 

15.0 

1.2 

3 

2 

Kirov 

MBS 

0 

38.2 

13 

46 

MBS 

0 

32.2 

5 

LMB-58 

1.79 

29.3 

15 

LMB-P 

15.0 

1.2 

13 

3 

Esheri 

MBS21 

2.91 

36.2 

6 

17 

MBS 

0 

26.3 

1 

1.92 

25.7 

10 

4 

Podolsk 

MBS 

0 

32.2 

4 

16 

LMB-58 

2.0 

28.4 

LMB-58 

1.8 

29.9 

2 

5 

Oussouriysk 

LMB-58 

2.2 

15 

38 

MBL 

10.6 

16.4 

11 

MBB 

o 

16.6 

12 

6 

Petropavlovsk- 

LMB-58 

22.7 

jgamm 

27 

IMB 

5.9 

tSmm 

7 

Yuzhno-Sakhalinsk 

MBS 

0 

32.2 

mmm 

32 

LMB-58 

9.0 

28.5 

imam 

8 

Zhvanets 

0 

32.2 

10 

24 

0 

32.2 

7 

LMB-58 

1.7 

23.8 

3 

LMB-58 

1.3 

32.3 

4 

9 

Mayly-Say 

LMB-58 

4.2 

waam 

31 

LMB-P 

1.2 

amm 

10 

Semipalatinsk 

MBS 

0 

32.21 

5 

27 

LMB-58 

1.6 

25.3 

11 

MB-10 

1.8 

17.2 

11 

11 

E35HB3SBHH 

IMB 

12.2 

5.9 

6 

6 

12 

MBS 

0 

30.8 

2 

2 

13 

Kooldoor 

EDMB-I 

0 

180 

3 

6 

EDMB-II 

0 

3 

14 

Kazakhstan 

MBS 

0 

28.4 

5 

5 

15 

Kap.  Yar 

MBS 

0 

33.0 

5 

5 

16 

MBS 

0 

29.4 

17 

17 

17 

MBS 

0 

34.2 

3 

3 

0  -  total  number  of  records,  prepared  to  digitization  and  processing. 
21  -  with  noise-reducing  arrangement. 
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Appendix  B 

Waveform  Parameters  and  IDG  Measurements 
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-  time; 

-  frequency,  [Hz]; 

)  =  2n  f  -  angular  frequency,  [rad/s]; 
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